Vortex Fluid Relaxation Model for Torsional Oscillation Responses of Solid 4 He 
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A phenomenological model is developed to explain new sets of detailed torsional oscillator data 
for hep 4 He. The model is based on Anderson's idea of a vortex fluid(vortex tangle) in solid 4 He. 
Utilizing a well-studied treatment of dynamics of quantized vortices we describe how the "local 
superfluid component" is involved in rotation (torsion oscillations) via a polarized vortices tangle. 
The polarization in the tangle appears both due to alignment of the remnant or thermal vortices 
and due to penetration of additional vortices into the volume. Both are supposed to occur in a 
relaxation manner and the inverse full relaxation time r _1 is the sum of them. One of them is 
found to change linearly with respect to the rim velocity V ac - The developed approach explains the 
behavior of both NLRS and AQ _1 seen in the experiment. We reproduce not only the unique V ac 
dependence, but also obtain new information about the vortices tangle, namely a divergence in r at 
extrapolated T~30 mK. 
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After the first report on "non-classical rotational in- 
ertia" (NCRI) in solid 4 He samples fl, the confirmation 
came from several torsional oscillation(TO) experiments, 
including some by the present authors 0. This finding 
had been discussed in connection to the iVCi?/Q of a su- 
persolid as originally proposed by LeggettQ. The mea- 
sured drop of the period is expected to appear due to the 
reduction of momentum of inertia, which originates in the 
appearance of a superfluid component which does not fol- 
low the rotation of the sample cell wall. It was, however, 
at problematically high T where the "transition" was re- 
ported. The following findings of the real onset tem- 
perature T [j| and possible vortex fluid (VF) state below 
T [1,0], would overcome the too high T c for BEC. And 
recent our workQ demonstrated a transition from the VF 
state to the supersolid state(SS) in solid 4 He below about 
75 mK. In the present work we discuss the vortex dy- 
namics in the VF state, including the supercooled condi- 
tion or measurements under " equilibrium conditions" 0] • 
On the other hand, there have been various measure- 
ments by now showing that this phenomenon depends 
on various parameters like pressure, impurity, sample 
quality H, " orientation" etc. Dependence of the phe- 
nomenon on the excitation amplitude of the TO (or on 
the rim velocity Vac) gives rise to special interest as pre- 
viously discussed [5j. Fig.l shows our TO results for the 
V ac dependence of both NCRI (namely relative change 
of period AP/P) and dissipation AQ^ 1 (inverse quality 
factor) of the 49 bar sample, the former(lower column) 
was given in our experimental work[5j and discussed as 
evidence for the VF state, pointing out the logV ac linear 
dependence, originally proposed by Anderson[6|. In addi- 
tion, the unique feature that signal decreases when V ac is 
increased, is argued to be evidence for thermally excited 
vortices in the VF state. Let us discuss some prop- 
erties of the observed phenomena^ depicted in Fig.l. 
(T).It is easy to see that the V ac dependence of the pe- 
riod drop disappears for some characteristic velocities 
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FIG. 1: New Data set of TO Responses throughout the vortex 
fluid state[||], including supercooled condition. Upper column 
indicate energy dissipationAQ^ 1 and the lower column shows 
nonlinear rotational susceptibility, N LRS=AP / APi oa d for 49 
bar hep 4 He at different T's as functions of V^ c - 



V ac < 10 - 30 fim/s. ©.The drop of period (NLRS) de- 
creases as the applied V ac increases. This is a sign that 
the superfluid part is being gradually involved in rota- 
tion. (3). For steady rotation with velocities exceeding a 
characteristic velocity this effect vanishes, the sample ro- 
tates as a whole. (4). The characteristic value of the ratio 
^-/AQ^ 1 is T and the pressure dependent quantity of 
order of unity. ©.At high T, the dissipation AQ _1 is 
a monotonically decreasing function of V ac , whereas at 
low T there is an obvious maximum. ©.One more fea- 
ture among reported results is the frequency dependence 
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of both P and AQ _1 [l(|. In the literature there was 
a speculation that this behavior can be associated with 
quantized vortices. For instance, in ll| it is pointed 
out that velocity V ac rj 10 \imj ' s coincides with the ve- 
locity created by a single circulation around the sample. 
Prokof'evjl] pointed out that "To understand why NLRS 
decreases with V ac , one has to consider the non- linear re- 
sponse of vortex loops and pinned vortex lines to the 
flow" . Huse et al. [13] developed a simple phenomenolog- 
ical model which introduces dissipative relative motion 
of two components realized via "phase slips" of quan- 
tized vortices. P.W. AndersonQ describes the scenario 
of a set of chaotic vortices (vortex fluid or vortex tangle) 
which under the torsional oscillation(TO) behaves like 
vortex-anivortex pairs in the Kosterlitz-Thouless model, 
and free (unbalanced) vortices bring the superfluid part 
into rotation. The role of vortices in rotation and tor- 
sional oscillation of solid helium had been discussed also 
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In the following, we propose a phenomenological model 
describing the behavior of the torsional oscillations in the 
presence of a vortex tangle. In a vortex free sample or in 
the case of an absolutely isotropic vortex tangle, the su- 
perfluid fraction does not participate in rotation or tor- 
sional oscillations. Therefore the momentum of inertia 
Ifull acquires a deficit I$f = p s VR 2 /2 where p s is the 
'local' superfluid density for the VF state, and V is the 
volume of the sample. Angular momentum of the super- 
fluid fraction appears only due to the presence of either 
aligned vortices (vortex array) or due to the polarized 
vortex tangle having nonzero total average polarization 
P =£ (s' z (£)) along the applied angular velocity fl (axis 
z, the magnitude of £1 = V ac /R, where R is radius of the 
sample). Here C is the vortex line density (total length 
per unit volume), s(£) is the vector line position as a func- 
tion of label variable £, s'(£) is the tangent vector. In the 
steady case there is a strictly fixed relation between the 
total polarization £(s'(£)) and applied angular velocity 



f2 = kP/2 = kC (s'(£)) /2. 



(1) 



Here k is the quantum of circulation. In a case when vor- 
tex filaments form an array the quantity C coincides with 
2D density n and fl} transforms to the usual Feynman's 
rule. Angular momentum of the superfluid part can be 
written as 'M.sf =Isf^ =Isf kP/2. 

The situation drastically changes in a nonstationary 
(transient or oscillating) case. The total polarization 
P(i) changes in time owing to both the vortex line den- 
sity C(t) and the mean local polarization (s'(t)} change 
in time according to their own, relaxation-like dynam- 
ics. Therefore the angular momentum of the superfluid 
part is not M.sf=Isf^ anymore. Because of relaxation 
processes there is retardation between fi(t) and Msj?(t), 
and the connection between them is nonlocal in time and 
Mgir(i) is some functional of time dependent angular 



velocity fi(t). There are two possible mechanisms for 
relaxation-like polarization of the vortex fluid. The first 
is an alignment of elements of the vortex lines due to 
interaction with the normal component (See [H| for de- 
tailed explanations). This interaction (mutual friction) 
is proportional to the local normal velocity , which in 
turn is proportional to the rim velocity V oc . Thus it 
is natural to suppose that polarization P of the vortex 
tangle due to alignment of filaments along fi(t) occurs 
with typical inverse time T ] _1 (V ac ) which is proportional 
to the rim velocity V ac . Let us illustrate the above with 
consideration performed in [IB} . In the presence of mu- 
tual friction there is a torque acting on the line and the 
angle 4> between axis z and the line element changes ac- 
cording to the equation dtft/dt = a(V ac /R) sm<j) (a is 
the friction coefficient, dependent, in general, on T and 
pressure p). Except for a short transient, the solution 
to this equation can be described as a pure exponential 
~ exp(— £/ti(V oc )), with the velocity dependent inverse 
time r-f (Vac) ~ aV ac /R. Thus, we conclude that dur- 
ing time-varying rotation or torsional oscillation vortex 
filaments tend to align along the angular velocity direc- 
tion. However, there can be not enough pre-existing vor- 
tex lines in the tangle to involve all the superfluid part 
into the rotation to satisfy the relation (JTJ , or on the con- 
trary the initial vortex tangle can be excessively dense. 
In this case deficient (extra)) vortices should penetrate 
into (leave from) the bulk of the sample. This penetra- 
tion occurs in a diffusion- like manner [161] and leads to the 
relaxation-like saturation of the vortex line density C(t) 
We assume that this saturation occurs in an exponential 



manner with some characteristic inverse time r 
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Due to linearity of the diffusion process we suppose that 
coefficient /3 is velocity independent, but can be a func- 
tion of T and p. Combining both mechanisms we assume 
that the whole polarization of the vortex fluid occurs 
in the relaxation manner with pure exponential behavior 
ip(t'/r) ~ exp(i'/r), and the inverse time t _1 of relax- 
ation is just the sum of r£" (V ac ) and t% , 

t- 1 =a(T)V ac /R + /3(T) (2) 

In the presence of relaxation the angular momentum 
M(t) of the superfluid part is related to the applied an- 
gular velocity fl(t) by the nonlocal relation, 
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Relation ([3]) implies that the angular momentum M(i) 
depends on the applied angular velocity H(t) taken in the 
all previous moments of time with the weight exp(— f/r). 
To clarify the physical meaning of constants a and b we 
consider the limiting cases of very small and very large 
frequencies. In case to — > the slowly changing func- 
tion H(t — t') can be considered as a constant and be 
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taken out of the integral, whereupon the rest of integral 
becomes unity and we have M w ^ = (a + 6)0. But at 
the same time, both components participate in the solid 
body rotation, thus (a + b) = If u u- In the opposite case 
of very large frequencies, lj—kx), the integral from rapidly 
oscillating functions fl(t — t') vanishes, so M^^oo = aVt. 
Since under these conditions the superfluid component 
does not participate in the motion at all, we conclude 
that the constant a is nothing but the full moment of in- 
ertia J/v of the sample without the superfluid part (which 
includes momentum of inertia of the empty cell I empty )■ 
Thus, the quantity b is moment of inertia Igp of the su- 
perfluid part. Substituting j3| with a — In and b = Isf 
into the equation of motion of the TO, we get 
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dt 
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i N n{t) + i SF / n(t -*>(-)— 
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+ k0 = O. (4) 



Here 9{t) is the angle of rotation of the oscillator, k is 
the spring constant. Relation^ is an integro-differential 
equation and, in general, not easy to solve. Because v?(^r) 
is a pure exponential function we can eliminate the inte- 
gral term. Omitting details we arrive at the case where 
equation ([4]) is reduced to an ordinary differential equa- 
tion of the third order, which has a solution in the form 
6(t) = 9q exp(iwi). The frequency to satisfies the relation 
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Thus, the frequency of the oscillation consists of three 
parts. The first one too = y/k/If u ii describes the os- 
cillation with full moment of the inertia If u ii as if all 
ingredients (empty cell, normal part, superfluid part) 
fully participate in motion. The second term is respon- 
sible for increase of the frequency because the superfluid 
component participates in the torsional oscillation only 
partly. The third term is the imaginary one. It describes 
the attenuation of the oscillation amplitude, i.e. it de- 
scribes the dissipation. The amplitude decreases (with 
time) as exp [— SJ(a>)i] , and the inverse quality factor is 
Q-x _ 2^{w) ^ Using the smallness of the Isf << I full 
we put u> = io full in the right hand side, yielding (index 
in LOf u ii is omitted) 
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Relations ©,((6]) are the final solution to the problem 
of the torsional oscillation when the superfluid compo- 
nent is involved in rotation via polarized vortex fluids, 
and polarization occurs in the relaxation-like manner. 
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FIG. 2: Parameters a(T), /3(T), and p 3 (T) obtained from the 
data of Fig.l with the use of analysis described in text. /3(T) 
goes to zero, or r to infinity at extrapolated T~ 30 mK. 



Being phenomenological, the approach developed does 
not allow determining some quantities entering the for- 
malism. Thus the parameters a(T,p) and (3(T,p) re- 
sponsible for the relaxation of the vortex tangle should 
be also obtained on the basis of the approach describing 
dynamics of quantized vortices, which is so far absent. 
Nevertheless comparison of our results with the experi- 
mental data allows us to explain a series of experimen- 
tal results and to get some quantitative information and 
insights. Let us analyze relations (0 and ©. From re- 
lations (O,© of our paper it follows that ^-/AQ^ 1 is 
equal to (1/2) (wr). It can take any value depending on 
the arrangement of the experiment. But usually obser- 
vations should be under conditions with lot on the or- 
der of unity. Therefore in many experiments the AP/P/ 
and AQ -1 are of the same order of magnitude, alth oug h 
sometimes they can be significantly different (see [12|). 
Dividing the first relation by the second one and taking 
the zero V ac limit in the relation A£y 'AQ -1 = (1/2) (cot) 
we get an expression for relaxation time (3{T) due to dif- 
fusion of vortices. Taking further the zero V ac limit for 
the period drop, and assuming that (3{T) abruptly van- 
ishes below the 'critical velocity' (which is equivalent to 
absence of vortices) , we find the superfluid momentum of 
inertia Isf, and, consequently superfluid density p s can 
be extracted from the graphs for AP/P. Knowing Isf 
(p s (T)), 0{T) and fitting the curves AP/P as functions 
of V ac it is possible to determine the inverse relaxation 
time due to aligning T 1 _1 (V ac ) ~ a(T)V ac / R and quan- 
tity a(T). Performing all procedures described above, 
we have all the necessary data, as shown in Fig. 2 where 
parameters ta(T), f3(T), and p 8 {T) are depicted. 

In Fig. 3 we show AQ^ 1 and ^p-=NLRS as functions 
of V ac , drawn using relations ([5]) and ^ and extracted 
experimental data. It can be seen that shapes of curves 
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FIG. 3: Energy dissipationAQ _1 and nonlinear rotational sus- 
ceptibility NLRS at different T as a function of V ac , obtained 
using relations (|5)l,(|6} with parameters taken from Fig. 2 



f_= lu/2-k dependence of NLRS and AQ" 1 observed in 
[lfj. Indeed, the significant dependence on lu appears 
when the inverse time r _1 of relaxation is comparable 
with uj, which can happen at higher T. In this range of 
parameters, the AP/P ([SJ is a monotonic function of u>. 
In summary the phenomenological model of relaxation 
processes of the VF state has been introduced. Unsteady 
rotation and torsional oscillation have been studied. De- 
pendence of both the NLRS and the AQ^ 1 on T, V ac 
and / have been studied. The results obtained may serve 
as a good qualitative description of the corresponding 
measurements in the VF state in solid 4 He. Combining 
theoretical predictions with experimental data it became 
possible to obtain some quantitative results. Actually re- 
cent experimental results [I3| can be well understood in 
terms of the present VF analysis, as an alternative to the 
interpretation in terms of superglass, by other authors. 
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and their response to the change of T correspond to the 
curves shown in Fig.l and Fig. 2. It is seen that in the 
limit V ac — > 0, or ui — > oo, or a(T) —> 0, NLRS reaches 
the maximum value. Physically it is clear, since under 
these conditions the superfluid part cannot participate 
in rotation at all. Other limits V ac — > oo, or uj — > 0, 
or a(T) —> oo correspond to the vanishing of the effect, 
which is also reasonable since under these conditions the 
superfluid part participates in the solid body rotation, 
and no effect appears. If relaxation due to diffusion (pen- 
etration) is weak, then the dependence of dissipation be- 
comes non-monotonic. Analysis shows that the critical 
value of t 2 _1 = (3{T) is equal to the frequency w. In fact, 
this can differ by some factor on the order of unity. One of 
the possible reasons for this difference is that we calculate 
using a purely exponential relaxation process, whereas in 
reality it can be better described by a more complicated 
dependence. The maximum value of dissipation AQ~^ ak 

should be at \ =rr and it should be reached at values of 
the rim velocity V QC = R{oj-(3{T))/a{T). This tendency 
is easily seen in Fig. 1, AQ~^ ak decreases with T and 
shifts in the direction of small V ac , then for some "criti- 
cal temperature" when r^ 1 = /3(T), the peak disappears 
entirely. It happens at T about 120 mK. Comparing 
with the experimental data one can conclude that the 
behavior described above indeed takes place for T above 
about 75 mK, but the agreement fails for lower T . It is 
remarkable that 75 mK was detected by authors of the 
present paper, as T c below which a hysteretic behavior 
takes place as a sign of a transition to a supersolid(SS) 
state(see Q). Relations ([5]) and © can also explain the 
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